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Abstract. Expansion for the free energy functionals of the Potts spin glass models 
with 3, 4 and 5 states up to the fourth order in Sqap around the replica symmetric 
solution (RS) is investigated using a special quadrupole-like representation. The 
temperature dependence of the IRSB order parameters is obtained in the vicinity of the 
point T = Tc where the RS solution becomes unstable. The crossover from continuous 
to jumpwise behavior with increasing of number of states is derived analytically. The 
' ^ I comparison is made of the free energy expansion for the Potts spin glass with that for 

tS ' other models. 

O 



> 

in 

m 



o 
o 



X 



PACS numbers: 64.60.De, 64.70.kj, 75.10.Nr 



Submitted to: J. Phys. A: Math. Gen. 



Potts spin glasses with 3, 4 o-nd 5 states 2 

The spin glass corresponding to the random variant of the Potts model [T] occured 
to be one of the first models where the scenario of the transition to the nonergodic 
state was discovered j2i |3j to be quite different from the standard one, that is of the 
Sherrington-Kirkpatrick (SK) model |4j. Nevertheless, it remains to be in the focus of 
investigations till now and a number of questions still are not answered (see, e.g., the 
recent papers Ref. |5i i6j that are also dealing with the mean field Potts glasses with few 
number of states). 

It is well known that in the case of SK model the so calUed full replica symmetry 
breaking (FRSB) [7] takes place at the instability point T^ of RS solution. All 
other models with two-particle interaction and with refiection symmetry behave in 
the same way [8]. The Potts glass models (with p > 2) are usually considered to 
be classical examples of models without refiection symmetry: FRSB does not work 
at Tc [21 [21 [9l [ini [m [12] and there is a region where the IRSB solution is stable 
(see, however, Ref. [5]). It is known (see, e.g., Ref.pJLj and references therein) from 
the numerical solution, that if the number of states p < 4, then the IRSB solution 
appears continuously at T^ - the bifurcation point of the RS solution. For p > 4 the 
IRSB solution appears jumpwise. The Potts model with p = 4 presents a boundary 
between two types of phase transitions. In this case the order parameter m defining the 
number of groups in IRSB construction achieves its maximum value m = 1 just at the 
point Tc = TinsB- It is worthwhile to follow in an analytic way, using the bifurcation 
theory, how the crossover from continuous to jumpwise behavior of the order parameter 
happens. This paper is devoted to this problem. To some extent we proceed in a same 
line as the authors of the paper Ref. [13] where the SK model was considered. 

It is worth to notice that the Potts symmetry leads to the cancellation of some terms 
in the free energy. This by-turn results in a trivial solution for the order parameters in 
the RS approximation. This is why the Potts models differ from other models without 
reflection symmetry, in particular, from the models similar to SK model, that is with 
the Hamiltonian 

H = -\Y.J^^U,U,, (1) 

with quenched interactions Jjj distributed with Gaussian probability, but without 
reflection symmetry: 

Tr \U^^^+^^] ^ (2) 

for some integer k. In this case the disorder smears out the phase transition: there is 
no trivial RS solution at high temperature (see, e.g. [T¥t [T3]). 

At the same time the Potts symmetry leads to the fact that the order parameter 
is one and the same for all of p states, that makes it related to two-particle generalized 
SK-type models and allows to show that the dependence of the free energy on the order 
parameters in the vicinity of the bifurcation point Tc is universal in both cases. 

We use a special representation of Potts models which originates from the analogy 
between three-state Potts spin glass and isotropic quadrupole glass (see below). In 
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our representation the number of integrations is less than in usual approach and so it 
is more convenient for actual calculations. Moreover, this representation demonstrates 
explicitly the appearence of cubic terms in the free energy expansion. These terms play 
a crucial role for the order-parameter behavior. 

The Potts glass model is defined by the Hamiltonian 



a — 2 ^ 'Jij^€ri,CTj- 

where p is the number of states and variables 0"i, aj can take the values 0, 1, 
Here Jij are random interactions distributed with Gaussian probability 
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where the factor A^ insures a sensible thermodynamic limit. 

The 3-state Potts spin glass was described earlier in our papers (see, e.g., 
Ref . |il6l W7\ [T2] ) using the representation in terms of the operators of quadrupole 
momenta: Q = fSJ^^ ~ ^)' ^ ~ ^ ("^^^ ~ '^y/i ^ = ^i Jz = ^i 0^ ~1- The operators 
Q and V commute and in the representation where both of them are diagonal have the 
following form 
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Throughout this paper we shall use the short-cut notations Qk 
Vk = (0, \/3, -\/3). 
The equation 

[QkQi + VkVi] + 2 = Q5u,i. 



-2,1,1) and 



(5) 



provides the equivalence of the Potts model with p = 3 to the model with the 
Hamiltonian 

h = -\y.[J^3Q^Q^ + g^,v,v,], (6) 

(with Jij = Gij). The quenched interactions are distributed with Gaussian probability. 

Let us note that this representation is very useful for actual calculations and allowed 

us to obtain, in particular, the number of metastable states at zero temperature [IT] 

and to determine the low-temperature boundary of the stability of the IRSB solution 

Using the replica method we obtain the averaged free energy of the system in the 
form 



F/NkT = lim - 

n^O n 



In Tr exp 



-2t^n + 



E {^^')' 



aj^P 



J2 (t^ {Q"Q^ + v'^v^ 



(7) 
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where g"^ is the glass order parameter, t = J/kT. Here indices a and /3 label replicas. 
The Potts symmetry permits to describe the glass state with one glass order parameter 
for both kinds of operators. 

Using the standard procedure we perform the first stage of the replica symmetry 
breaking (IRSB) according to Parisi {n replicas are divided into n/m groups with m 
replicas in each group) and obtain the free energy in the form (with g"^ = ri if a and 
(3 are from different groups and g"^ = ri + v ii a and (3 belong to the same group). 



^iRSB = -NkT X 
t^ I2 + {I - m)^ - 2vA + 

-\n fdz'' [ds^'l^r]. (8) 



m 



Here 



^ = exp (-2^1) + exp (^i) [exp (^2) + exp (-^2)] , (9) 

e^ = zt^, 62 = stVSv, (10) 




^" '...)expl--)= Idz^i...). (11; 



^27r 

The Potts symmetry provides also the trivial solution q^is = and ri = (see e.g. 
Ref [H]). Let us note that in general case of the Hamiltonian ([6]) with Jij 7^ Gij one of 
the glass order parameter can not be zero. This fact can be easily seen by analyzing the 
high-temperature expansions for the RS equations for the order parameters [18]. Let 
us note that the Eq.(l8]) can be obtained from the appropriate equation of Ref. \TT\ by 
the integration over suitable linear combinations of variables. 

Expanding the free energy (IHl) near Tc up to t;^ and solving the equations that define 
the extremum conditions for the free energy one can obtain [16j : 

v = 8t- 84r^ m = r, AF = 16r^ (12) 

where AF = Firsb - F^s and r = t -tc- 

Let us consider now the case p = 4. In this case one can use the following operators 
Q' ~ 3Jz^ — J (J + 1), V ^ Jz for J = 3/2 and the third operator P' defined as to 
be orthogonal to Q' and V. Namely, 2Q'^ = (1,-1,-1,1), 2V5V^ = (-3,1,-1,3), 
2v^P^ = (1,3, -3,-1). 

It is easy to see that the following equations hold: 

1 12 12 

^ 4' 4 5^' 4 5^' 

V'P' = -—Q'; Q'V' = -—P' + -Q'- 
10^' ^ 10 5^ ' 

Q'p' = Ip' _ Aq' 

^ 5 10^' 
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so that 

Q'uQ'i + ViVl + PlPl + \ = 5,,. (13) 

There is no reflection symmetry and we have 

TrP'V = --; TrV'2g' = -; Ti P'V'Q' = - — , (14) 

that leads to the appearance of cubic terms in the free energy expansior|j|. 

Now, using the rephca method one can write the averaged free energy for the 
Hamiltonian 

h = -\y: hq\q', + y^i + Kp% (15) 

in the form analogous to ([7]) and perform the IRSB. The extremum conditions for this 
IRSB free energy have the form of three equations that look different only at the first 
glance. In fact they are the equations from Ref.[lT|. integrated over different variables 
and so have a common solution v = {{Q'aQ'p)) = {(^L^p)) = {{PL^'p))- Expanding the 
IRSB free energy in the neighborhood of tc up to f ^ we obtain (in terms of standard 
variables (as, e.g., in [TT])): 

AF/Nkt = — (-1 + m)t^v^ x 
8 

X [-6 + 6t^ + 4(m - l)t^v + (Sm^ - 3m - 7)t%% (16) 

Here t = tc + t, m = mo + 6m and r, v and 6m, are small. The extremum conditions 
relative to v and m give the system of equations: 

(-1 + m)t^v X 

X [-3 + 3t^ + 3(m - l)tS + (3m2 + 3m - 7)^^;^] = 0; 

-t\2[6 - 6t^ - 8(m - l)t% + (4 + 12m - Qm^)fv'^] = 0. (17) 

8 

To determine the relative orders of the small values r, v and 6m,, one has to use the 

so-called bifurcation equation (see, e.g. Ref.|T9]) for the system (IT7|) . One can also 

check directly that the equations (IT71) become incompatible if v and 6m, are supposed 

to be series in integer powers of r. If 6m, = mx^fr and v = v^^ (as it follows from the 

bifurcation equation) then we obtain tc = 1, "^o = Ij and for m\ and v\ the following 

system of equations can be obtained from ([1] 



6 + 3mit'i — lv\ = 0; 

12 + Smivi - 7vl = 0. 

I It is worth to notice that our representation in terms of three matrices is not unique. For example, 
one can use the set: Ak = (1,-1,— 1,1), Bk = (—1,1,-1,1), Ck = (1,1,-1,-1). We use, however, 
the quadrupole-hke operators in order to retain the physical meaning of generalized anisotropic models 
and to preserve the analogy between even-p and odd-p cases. 
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The system has two solutions that differ only by the sign. Taking into account that 

AF/NkT = -mlvfr^/^, 
and that on the physical branch m < f , we choose 



f2r/35; m=l- ^21t/5. 

So, in the case of Potts spin glass with 4 states m = 1 at the point of bifurcation tc, 
that is tiRSB = tc- 

Let us consider now the Potts glass with p = 5 in an analogous way. We shall 
use the following diagonal operators: P^ = (—4,1,1,1,1), P^ = (0,-3,1,-1,3), 
pf ) = v/5(0, 1, -1, -1, 1), Pf ^ = (0, 1, 3, -3, -1). The equivalence to the Potts model 
follows from the relation: 

Pt'^Pi'^ + Pi'^PF^ + Pi'^PP + Pt^P!"^ + 4 = 204,z. 
The nonzero cubic terms are 

p(l)^ = _12; p(2)2p(l) ^ 4. p(3)2p(l) ^ 4. 



p(4)^p(l) ^ 4. p(2)-p(3) ^ ^g/^g. 
p(4)2p(3) ^ _16/y5. p(l)p{3)p(4) ^ _12/y5 



Now the IRSB free energy expansion near tc has the following form (in standard 
variables) : 

AF/NkT = - — (-1 + m)tV x 

X [-12 + 12t2 + 4(2m - l)tS + (Gm^ + 6m - 35)t''t;2]. (18) 

The extremum conditions are 

(-1 +m)t^v X 

3^ ' 

X [-6 + 6t^ + 3(2m - l)t^v + (Gm^ + 6m - 35)tV] = 0, 

—t^v'^[l2 - 12t2 - 4(4m - 3)t% + (41 - 18m2)tV] = 0. (19) 

It is easy to obtain that now v and 5m can be presented as series in integer degrees 
of r: V = ViT + f2T^, m = mg + niiT. Here t;i, ^2, mo and mi satisfy the following 
equations: 

6 — 3t'i + imQVi = 0, 



and 



4 — f 1 + 2mof 1 = 

6t;2 - 12mi - 92 = 0, 
6f2 — 16mi — 43 = 0, 
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Figure 1. (Colour online) Schematic representation of the order parameter v in the 
vicinity of the bifm'cation point tc for p — 3,4,5. Solid lines for tti < 1, dotted lines 

for TTl > 1. 



SO that finally we obtain 



tr 



1; m 



3 49 

— I t: V 

2 4 ' 



-2r 



239 
~6~' 



and the transition from RS to IRSB can not take place at the bifurcation point. 

Here we would like to make a remark about the 5-state model without any pretend 
to be rigorous. Basing on the fact that the numerically obtained \\A\ value of ti^sB 
for m = 1 is very close to our tc, one can hope to obtain the jump to the IRSB 
solution already from the equation ( ITSl) and to estimate qualitatively the characteristics 
in question. In fact, if we expand (IT8|) in the vicinity of m = 1 and let its first and 
second derivatives relative to v be equal to zero, we obtain t\^gQ = 184/187 and the 
jump Av = 561/8464 (compare with the table in [11]). 

Qualitative behavior of the solutions for v near tc for p = 3, 4, 5 is presented on the 
figure [H 

Now let us consider the IRSB free energy expansion near the bifurcation point tc 
from a more general point of view. It is interesting that the form of the series for AF 
in small deviations 6q°'^ from qjis up to the third order is one and the same for different 
models and coincides with the one that can be written [15] for the random generalized 
Hamiltonian ([1]) where U is arbitrary diagonal operator. The fact is that in the most 
general case (including ([T]) and (J3])) AF depends on the following (the only nonzero) 
sums: 



lim -^ {Sq"^f = -[r-{m- l)v]^ - m(l 



n^O n 



m]v 



a,/3 



hm - ^ 'dq^'^Sq"^ = [r - 



n-i>0 n 



m 



l)v] 



(20) 
(21) 



«,/3,(5 



lim — y^ {5q 

a,l3 



a/3\3 



[r - {m- l)v]^ + 
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?im{m — l)[r — (m — l)!"]!!^ + m{m — \){2m — l)v^, (22) 

lim - y 'dq'^^dq^^dq^'' = 2[r - (m - l)tf - 

3m{m — l)[r — (m — l)f]f^ — m^{m — l)v^, (23) 

lim - y ' idq'^^Ydq''^ = W - (m - l)vf - 

a,p,'y 

[r-{m- l)v]v^, (24) 

lim - y '6q"^6q"^6q"^ = 

lim - y '5q''^5q'^^5q^^ = -[r - {m - l)vf , (25) 

with r = ri — q^s- The prime on the sum means that only the superscripts belonging 
to the same 5q are necessarily different in Y^' . 

So, the deviation AF of the free energy FmsB from its replica symmetric part in 
the most general case is 



AF t^ 

NkT 4 


1- 


-ew 


( 

\ — ^ — {m — \)v\ 


v^m{l - m)| - 


-tL[r-{m- l)vf - 


t"" [C [r - (7 


71 - 


-l)v] 


+ 



D\r — (ra — \)v\ v^raim — 1) — B^v^m^^m — 1) + 
B4V^m{m-l){2m-l)} + ... (26) 

where t = t^ + r, and parameters W, L, C, D, B^, B4 are some combinations of operators 
averaged over the RS solution. For example, the coefficient L enters AF with the sum 

Let us note that in the case of zero RS solutions for the order parameters the 
expansion does not contain the terms where some indices occur only once. In the case 
of reflection symmetry there is no terms where some indices occur odd number of times. 
If there is no reflection symmetry, using the extremum conditions for the free energy 
for Hamiltonian (jS]) and taking into account the fact that L\t=t,. 7^ 0, the bifurcation 
condition gives|15] r — {m — l)v = + o(At)^, i.e. the condition that there is no linear 
term for the glass order parameters. In fact, there is no other linear term because 

1 - t^W] \t=t. = ARSrepl|t=t. = (27) 

at the bifurcation point. 

In the case of the Potts spin glass model the situation is quite different, although 
the reflection symmetry is absent. Now L = because it is zero RS solution that 
bifurcates so that single indices mean that L is composed of zero values of averaged 
TtQ,TtV.... Analogously the coefficients in front of ( 12^ and ( l25l) are also zero. So, we 
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have C = 2B^ — B4 and D = 3{—B^ + B^). Here S3 is the coefficient in front of ,1^^,, . 
and i?4 that in front of (122|) . Since -^|t=jo = 0, then the condition r — {m — l)v = is 
not fuffilled. (If J9 > 2 then ri = r = 0.[llj) 

In general case we obtain from the extrenium conditions for the free energy (126|) 
the following equations for the order parameters: 

ttmirri 



2m{m — l)vZAt 



l)v{3[-B4+ 



m( 



-B3 + 254)] v + 2D[r-{m- l)v]} 



(28) 



(2m - l)v^ZAt = t^v^ {(2m - 1) [-^4+ 
m(-53 + 254)] V + m{m - l)(-53 + 2^4)^; + 
D{2m- 1) [r - (m- l)v]} , 



where 



d 


"t2 




" 






1- 


-t'^w\ 


di 


_T - 





Hence we obtain from (l28|) - (l29|) : 

m{m — 1) [mBs — B4] = 0. 
In fact, in the case of the Potts spin glass with p 



(29) 
(30) 

(31) 
3 for the Hamiltonian (161) at 



tc we have: 
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B, = -r 

D 



12 
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(TrQ^Sy + {TrVy3 
TrQy3y + 3(TrQVy3 



(32) 



In this case in the neighborhood of the bifurcation point one can get from the 
equations (p8|) - (l29|) that f ~ r, and from the equation (13T|) m = |. 

Analogously, in the case of 4 states (for the Hamiltonian (IT5|) we obtain: D = 0, 
i?3 = -B4 = Y28 ^'-^ ^^^^ rri = 1. The equation (l28l) is the identity now. The r.h.s. of (l29ll 
becomes zero and one has to take into account higher order terms in the expansion of 
FiRSB- This leads to t; ~ y/r. For p = 5 we have D = 64, B3 = i|^,i?4 = 64, so that 
m = I and v ~ — r. 

It is worth noticing that from the free-energy expansion (|26l) it follows directly 
that the stability of the RS state is determined by the sign of A^sreph just as in the 
case of SK model |20| . To prove this fact let us consider small deviations from the RS 
solution: q°''^ = qns + ^Q""^ at arbitrary temperature and let us perform IRSB. Using 
the notations introduced before, (that is Sq'^^ = r if a and /3 are from different groups 
and 5g"^ = r + f if a and /3 belong to the same group) we can write Firsb — Frs up to 
the second order in 5q in the form (126]) with arbitrary t. Since [1 — t^W] = ARSrcpb it is 
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easy to see that for m < 1 and L > the RS state is stable (AF < 0) for the values of 
temperature such that A^srcpi > 0. The similar result is valid for all subsequent stages 
of RSB because of the Parisi rule 

\im -y' (5q''^Y <0 

To conclude, we introduced the representation for the Potts glass model based on 
the operators of the quadrupole momenta. With the help of this representation the 
expansions for the free energy functionals of the Potts spin glass models with 3, 4 and 
5 states up to the fourth order in Sqa/s around the replica symmetric solution were 
obtained. The temperature dependence of the IRSB order parameters in the vicinity of 
the point T = T^ where the RS solution becomes unstable was derived. The crossover 
from continuous to jumpwise behavior with the growing of the number of states was 
traced analytically. The comparison was made of the free energy expansions for the 
Potts spin glass with that for other models and the similarity of such expansions was 
demonstrated. 
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